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^ ■ Abstract 

o ■ 

The coefficients that describe the anisotropy in the azimuthal distribution of particles 
are lower when the particles are recorded in a detector with finite granularity and measures 
only hits. This arises due to loss of information because of multiple hits in any channel. 
The magnitude of this loss of signal depends both on the occupancy and on the value of the 
coefficient. These correction factors are obtained for analysis methods differing in detail, and 
■ are found to be different. 

<N 
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1 Introduction 

in . 

Azimuthal anisotropy in particle emission in ultra-relativistic heavy ion collisions was 
proposed as an important probe of the dynamics of the system [1]. Subsequently, 
various methods have been proposed to obtain this anisotropy [2, 3, 4]. In the more 
commonly used method the azimuthal distributions are expanded in a Fourier series 
where the coefficients of expansion are the measures of different orders of anisotropy 
[2]. For small values of these coefficients, the first two terms describe an elliptic shape. 
The first order anisotropy v±, the directed flow, measures the shift of the centroid of 
the distribution and is the coefficient of the first term in the expansion. The second 
order anisotropy t>2, the elliptic flow, measures the difference between the major and 
minor axes of the elliptic shape of the azimuthal distribution and is the coefficient 
of the second term in the expansion. The elliptic flow, t>2, probes the early stages 
of expansion of the interacting system and has been measured by large number of 
experiments for different particle species in different kinematic domains for a variety 
of colliding systems and a range of center of mass energies [5, 6, 7, 8, 9, 10, 11, 12, 13]. 
These measurements have provided new perspectives on the observed mass dependence 
of the elliptic flow [14]. 

In detector sub systems where pt or energy is not measured, the anisotropy coeffi- 
cients are determined from the azimuthal distribution of the number of particles. Some 
detectors measure the distribution of hits 1 . If a cell is hit by more than one particle, 
information is lost because the cell is still registered as one hit. The (At/ x A<3>) size of 
each cell in Silicon Pad Multiplicity Detector in WA98 experiment is about 0.07 x 2° 
[15], whereas the corresponding size in the silicon strip detector in NA50 experiment is 
about 0.014 x 10° [16]. For the same occupancy, both detectors will lose comparable 
number of particles by measuring hits. Therefore, the anisotropy coefficients describing 
distribution of hits are expected to be smaller than coefficients describing distribution 
of particles; v^ ts < v n and there is need to determine an appropriate correction factor 
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1 One activated cell is counted as Nhits = 1 irrespective of the number of tracks activating it. 



which will depend on the granularity of the detector. Methods to estimate the effect 
of finite granularity have been discussed in [17]. In the present work the effect of finite 
granularity on the standard methods of analysis is investigated by folding the detector 
geometry in the simulated data. 

In Section 2, an approximate expression for the ratio ^ — is obtained as a function 
of occupancy. Section 3 describes the simulation and the different methods of analysis. 
The Results are discussed in Section 4. 



2 Multiple Hits and Azimuthal Anisotropy 

It is possible to have an ideal setup for a simulation experiment for any conceivable 
granularity. However, an actual experiment has a finite number of detector cells N ce u s , 
which defines the coarseness of the granularity for a given acceptance of the detector. 
For a given average number of incident particles, < N part >, one can define the mean 
occupancy /io 

^ Npart ^* / -■ \ 

V0=— (1) 

^cell 

Using the Poisson distribution for probability of n particles incident on any cell, 
one can deduce the average number of hits as 

< Nms > = 1 ~ e~ M " = -x (2) 

< N part > A*0 ln(l-x) 

where x = < Nh its > /N ce u is the hit-occupancy and is experimentally measurable. 
Mean occupancy can also be written as /io = —ln(l — x). Since the total number of 
cells (N ce u) is the sum of the average number of occupied (< N occ >) and unoccupied 
(< N unocc >) cells, one can immediately obtain the expression for occupancy as in 
ref [13]. 

= Mi + ^° CC> J (3) 

^ unocc 

enabling its determination from experimentally measurable quantities. 

The anisotropy in the azimuthal distribution of the number of incident particles is 
written as 

N part {<t>) =< Npart > + 

cosn(0 — ip n )) (4) 

where ip n is the event plane angle. Using equations 1 and 2, one can write the 
azimuthal dependence of hits as 

iW0)ocl-e-^<« (5) 

where denotes the azimuthal dependence of the occupancy. Since the intrinsic 
occupancy of cells increases with the increase in the number of incident particles, the 
occupancy will have the same azimuthal dependence as the incident particles and can 
be written as fj,((f>) = /x (l + 1] 2w n cosn(0 — ip n )). Substituting this in equation 5 and 
expressing Nh its as a Fourier series with coefficients v^ ts enables a determination of 
the ratio v^ ts jv n . To the first order, this ratio can be approximated as 



where f(v) = fi v 2 for n = 1, and = for n = 2. 

The function f(v) contributes little for small values of occupancy and flow. 
The ratio can also be approximated as 

= ■ ln(l - x) (7) 

v„ X 

These results have been applied to the data recorded in the Silicon Pad Multiplicity 
Detector in the WA98 experiment [8]. The results from equation 6 and 7 are 
corroborated with results from simulations as described in the following. 



3 Simulation and Analysis 

For the present simulation experiments, the events have been generated with various 
values of charged particle multiplicity corresponding to different occupancies in the de- 
tector. Assuming a constant dN/dr/ and an exponential p^ distribution, the kinematic 
variables of each particle are generated with p T in the region to 6 GeV/c and r\ in 
an assumed region of acceptance of the detector. Typical ranges of rj chosen in the 
present work vary between 0.5 and 1.0. Azimuthal angle of each particle is assigned 
according to the probability distribution [18] 

r(0) = [1 + 2«i cos(0 - + 2v 2 cos 2(0 - ^ R )\ (8) 

where i/jr is randomly generated once for each event. Events are generated for different 
granularities in rj and $. A constant dN/dry distribution and cells of equal A77 intervals 
give a uniform intrinsic occupation probability for each cell. Detector geometry, flow 
and occupancy are varied for a systematic study. 

In the present work, the granularities that are chosen are fairly arbitrary but com- 
mensurate with the coarseness of certain detectors [15, 16]. More specifically, simula- 
tions are performed for Ar] = 0.07, A$ = 2°; Ar] = 0.014, A$ = 10°; Ar] = 0.00875, 
A$ = 10°. The results are based on an analysis of 10 6 events in each case. For low 
flow values (v n = 0.02), the number of generated events is 2.5 • 10 6 . 



3.1 Different Data-Sets 

The number of particles simulated and the number of cells activated (Nh its ) are known 
for each event. It is assumed that one incident particle does not activate more than 
one cell. 

The following information is stored as three different data-sets from the simulated 
events. 

(a) the number of particles and the azimuthal angle of each particle. This corre- 
sponds to a measurement in an ideal detector of infinite granularity. Anisotropy 
coefficients measured thus are labeled as v^ eal . 

(b) the number of particles and the azimuthal angle of each hit cell. This corresponds 
to the case when azimuthal angle of each particle is known to an accuracy de- 
termined by the azimuthal size of the cell and the number of particles can be 
determined using the pulse height information. This is equivalent to randomly 
adding (or subtracting) 8(f) ( < 4p) to each where A$ is the azimuthal size of 



each cell. The anisotropy coefficients obtained using these are called and can 
be written as [19] 

rh ■ A<J> 
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(c) the number of hits and azimuthal angle of each hit cell. This is the information 
recorded by the detectors that produce only a binary (hit/no-hit) signal for each 
cell and the corresponding anisotropy coefficients are called v^ ts . 

3.2 Methods of Analysis 

Fourier coefficients of n th order can be determined from the azimuthal distribution of 
the particles with respect to the event plane angle of order to, provided n is an integral 
multiple of to, by fitting to the following equations [2] 

dN 

d{(f) _ y j oc 1 + E^ =1 2< m cos nm(0-<) (10) 
The event plane angle is given by 

C = -ftan- 1 p^f ] (ii) 

to \ hwi cos mcpi J 

where the summation is over all particles % and the weights Wi are all set to 1. The 
average deviation of the estimated event plane from the true event plane due to multi- 
plicity fluctuations can be determined experimentally and is termed as the resolution 
correction factor (RCF). Experimentally, RCF is obtained using the sub-event method 
described in reference [2]. Here every event is divided into two sub-events of equal 
multiplicity and the event plane angle ip' m is determined for each sub-event. This en- 
ables determination of a parameter Xm directly from the experimental data using the 
fraction of events where the correlation of the planes of the sub-events is greater than 
tt/2 [2, 20]: 

N events (mW«-^\>n/2) = ej^_ 

Ntotal 2 

where N to tai denotes the total number of events, i/j'^, ip!^ are the observed event plane 
angles of the two sub-events (labeled a and b) and the numerator on the left denotes 
the number of events having the angle between sub-events greater than 7r/2m. The 
parameter Xm is used to determine RCF„ m = {cos(nm(ip' m — ip 1 ^)))-, where ip^ l ue is 
the true direction of the event plane, and the average is over all events. The RCF can 
be determined from Xm by the following relation in reference [2]. 



(cos(nm(^ - ^ ne ))> = ^|x m exp(- X ^/4) • 

V0d/4) + 7^(^/4)] (13) 

where l v are the modified Bessel functions of order v. The RCF can also be obtained 
by obtaining (cos n(ip'^ — ip'n)), where ip'^ b are the event plane angles of the two sub- 
events. 



In the present work, Fourier coefficients v' nm are extracted for the case with event 
plane order equal to the order of the extracted Fourier coefficient, i.e. v' nm = v' nn and 
is denoted here by v' n . The values of v n have been obtained by the following methods: 

Method 1 : In this method, the sub-events for each event were formed by dividing 
the pseudorapidity range into two such that each sub-event has equal number 
of particles(hits) 2 . Then v'® = (cosn(0" — V'n)) an d v'£ = (cosn(<^ — VC)) are 
determined where 0" represent the azimuthal angles of particles in sub-event a 
and il>% is the event plane angle determined using particles in sub-event b. The 
averages are computed over all particles over all events. In the absence of non 
flow correlations 



v n = 



n n (14) 



^(cosn(^-^)> 



It is also possible to obtain v' n by fitting equation 10 to the <pi — ip' n distribution. 
This distribution is a sum of the the distributions 0" — ip'^ and (j>\ — ip'^. This 
yields 



(15) 



The denominator in both cases above is the event plane resolution correction 
factor when the event plane is determined for the sub-event with half of the 
complete event multiplicity and is approximately lower than the full event RCF 
by a factor \/2. The values of v n obtained this way are termed as v^ eom . 

Method 2 : In this method, in each event, the sub-events were formed by randomly 
selecting one half of the particles. The v' n values are also extracted by fitting 
equation 10 to the <p — ip' n distribution, where ip' n is obtained by excluding the 
particle (or hit) being entered in the distribution 3 . The v n values are obtained 
from 

^WF. (16) 

RCF n is the resolution correction factor for the full event plane and is obtained 
by the correlation between randomly divided sub-events of equal multiplicity and 
using equations 12 and 13. The values of v n obtained this way are termed as 

r) .rand 



4 Results and Discussion 

A relation for the ratio of the anisotropy measured using hits to the actual anisotropy 
for different values of occupancy was obtained in section 2. These values are corrob- 
orated using simulations and the results are discussed in this section. The different 
methods of analysing the data discussed in section 3 are applied on simulated data to 
study the effect of finite granularity on the v n values. The simulated data are analysed 



2 The two sub-events were also formed by assigning particles(hits) to each from alternate segments of the 
imuthally segmented detector. 

3 This avoids autocorrelations but also introduces a negative correlation, effectively decreasing the values of 



for all the detector geometries described above. The results are discussed for the case 
Arj = 0.07, A$ = 2°. The conclusions remain the same for the other geometries. 
For all simulations, ^hits/^part is obtained for different values of x and is found to be 
consistent with results from equation 2. 

4.1 v n Using Known Event Plane: Actual Dilution 

The dilution in the anisotropy coefficients due to finite granularity can be computed 
using the known event plane angle (ifj n ) in simulation. The quantity v n = (cos(n(0j — 
i/j n ))) is determined for the different data-sets described in Section 3.1 and yields the 
finite granularity effect on the anisotropies in the distribution. 

1. vlf eal reproduces the input flow, as expected naively. 

2. The dilution due to coarse information about the particle angle can be judged by 
plotting the ratio v 1 ^ jv % ^ eal . The result for the granularity Arj x A$ = 0.07 x 2° 
has been plotted in Fig. 1 for two different values of initial anisotropy. The 
correction factor due to coarse information of particle angle is very close to 1 for 
such a small azimuthal size of the cells. The results have been corroborated using 
simulations for values of A$ up to 30°, and agree with results from equation 9. 

3. The anisotropy for the hits, v^ tts is diluted both due to coarse information of 
particle angle and loss of particles because of the multiple hits. The resultant 
loss is best seen by plotting the ratio v^ ts /v % ^ eal as a function of hit-occupancy 
for the simulated data. The results are shown in Fig. 1 for a A$ = 2° along 
with the estimates obtained using equations 6 and 7. For A$ = 2°, results in 
(ii) above show that the coarse information of particle angle has very little effect, 
and the dilution in v^ ts is primarily due to loss of particles because of multiple 
hits. Simulation results corroborate the analytical expression that include a weak 
dependence on anisotropy v. 

4.2 v n Using Reconstructed Event Plane: Observed Dilution 

In this section we investigate the results on dilution of anisotropy when the event 
plane and its resolution is determined from the data. Both the methods listed in 
Section 3 require a determination of (i) the uncorrected v' n and (ii) the corresponding 
(sub) event plane resolution. The quantitative effect of finite granularity on each of 
these quantities is different, and hence the measured values of v n are different from 
the initial values. The results of a systematic investigation are shown in figure 2 for 
varying hit-occupancies. The results from method 1 of section 3 are shown in the left 
column and those from method 2 are shown in the right column. For both methods, 
the results from data-sets (b) and (c) are scaled by the corresponding values obtained 
using data-set (a). The open circles show the results obtained using the data-set (b) 
for all charged particles and the mean angles of the cell positions. The filled circles 
show the results obtained using the data-set (c) for the hits and the corresponding 
angles. For both methods, analysis of data-set (a) reproduces the initial anisotropy, 
validating the methodology. 



4.2.1 Division into sub-events based on geometry: vf^ om 

The anisotropy coefficients v^ eom determined using equations 14 and 15 yield identical 
results. 

The event plane resolution correction factor and the uncorrected values of v[ are 
seen to decrease by different factors due to finite granularity effect, resulting in a 
reduced value of vf eom . 

The results show that the value of anisotropy measured using a detector with 30% 
mean hit-occupancy is to be corrected by a factor of ~ 1.2 to obtain the actual value of 
anisotropy (i>i), clearly a significant effect. The results for the second order anisotropy 
are similar, with small quantitative difference as seen from equation 6. 

The analysis described as Method 1 was repeated for the case when the two sub- 
events were formed by assigning particles(hits) to each from alternate segments of the 
azimuthally segmented detector. The results remain the same. 

4.2.2 Random division into sub-events: v™ nd 

The right column of Figure 2 shows the results for method 2 when the events are 
divided randomly into two sub-events, and the projection of particles/hits is taken 
on the event plane angle of the full event (after removing autocorrelations). This 
method works for the simulated data corresponding to data-sets (a) and (b) described 
in Section 3.1, and the values of data-set (b) scaled to the corresponding values from 
data-set (a) are shown as open symbols in the right column of Fig. 2. The values 
obtained using data-set (c) are much lower than the values obtained by analysing 
the same data using Method I. The systematically lower values of the quantities for 
Method 2 arise due to multiple hits, due to removing auto correlations and due to 
random division into sub-events. The combined effect results in much lower values 
of v\ and . For data-set (b), the method of removing autocorrelation removes only one 
particle while the other neighbouring particles are used in determining the event plane. 
For data-set (c) comprising of hits, one detector cell is removed from the data that 
determines the event plane, effectively removing all particles within the azimuthal size 
of that cell, introducing a negative correlation, resulting in much lower values of v[. 
This holds true for all values of occupancy. 

The decrease in the values of event plane resolution can be understood as follows: 
For an azimuthal distribution given by equation 8, the particle density is maximum 
along the direction of the reaction plane. On an event by event basis, the maximum 
loss of particles due to multiple hits will be along this direction. Consider there are 
N corr correlated particles in a region 5(f) about the reaction plane, where N corr / N tota i 
is greater than 50/27T. When such an event is divided into two equal multiplicity 
sub-events, the correlation between the two sub-events will be maximum if N corr /2 
particles go into each sub-event. Though this is true on the average, on an event by 
event basis, only a certain number out of N corr fall into one sub-event. The correlation 
between sub-events for these events is less than the corresponding situation described 
above. The correlation will be weakest if all of these particles fall into one sub event. In 
such a situation, the v™ nd values are likely to be much lower than the v ldeal . However, 
the probability of the random division into sub-events leading to this situation is 
[\j2) Ncorr ~ 1 and is small. 

The situation remains the same when hits are recorded instead of particles, and 
N corr is replaced by N corr ( hits ), and N tota i by N hits . The probability of a random division 
with all N hits going into a sub event is (\/2) Ncorr{ - hit! >)~ 1 . This probability is clearly 



greater than the corresponding case where there is no loss of particles due to multiple 
hits. When this happens, the two sub-events show little correlation, causing both v' n 
and RCF n to decrease. The decrease in both quantities only partially compensate each 
other and the resultant values are observed to be lower, as shown in the third panel 
on the right column of figure 2. 

The resultant effect is to produce a reduced value of v™ nd , which would need a 
relatively larger correction factor to obtain the original value. Repeating the simula- 
tion for different values of anisotropy shows that the correction factor increases with 
decreasing values of anisotropy. Below a certain value of anisotropy, v thres , the values 
of anisotropy extracted using this method yields v[ and consistent with zero, putting a 
limit on the sensitivity of detecting anisotropy in the data. The values of v thres de- 
pend both on the granularity of the detector and on the order of the anisotropy being 
determined. 



5 Conclusions 

The anisotropy in the distribution of hits is shown to be lower than the anisotropy in 
the distribution of particles. Loss of particle information due to multiple hits (or two 
track resolution) contributes significantly to the dilution of the observed anisotropy 
values. The correction factor for the dilution has been obtained and confirmed by 
simulation experiments for different detector geometries. Different methods of event 
subdivision yield same results for an ideal detector. While it is known that the event 
subdivision obtained randomly does not work in the presence of non-flow correlations, 
the limitation of this method is shown here for the case where there are no non-flow 
correlations. The anisotropy values v ™ nd obtained using the distribution of hits are 
much lower than the corresponding values of v^ eom and need a correction factor which 
is larger than the one obtained in Equation 6. The analysis and the discussion in 
the present work are suited for the case where the event plane is being determined 
from the same set of particles. The event plane resolution correction factor corrects 
only for fluctuations arising due to finite multiplicity. If the anisotropy parameters are 
determined with respect to an event plane determined from another set of particles 
measured using a detector with different granularity, then the correction factors need 
to be determined differently and their determination is outside the scope of the present 
work. 

Acknowledgement 

The financial support from the Department of Science and Technology and the Depart- 
ment of Atomic Energy of the Government of India is gratefully acknowledged. 

References 

[1] J.Y. Ollitrault, Phys. Rev. D46 229 (1992). 

[2] A.M. Poskanzer and S.A. Voloshin, Phys. Rev. C58 1671 (1998). 

[3] N.Borghini,Phuong Mai Dinh and J.Y.Ollitrault Phys. Rev C64 054901 (2001) 

[4] R.S.Bhalerao, N.Borghini and J.Y.Ollitrault Nucl.Phys.A727 373 (2003). 

[5] E877 Collaboration, J.Barette et.al., Phys. Rev. Lett. 73 2532 (1994); Phys. Rev. Lett. 
70 2996 (1993); Phys. Rev. C55 1420 (1997); Phys. Rev. C56 3254 (1997). 
E802 Collaboration, L.Ahle et.al., Phys. Rev. C57 1416 (1998). 



[6] WA93 Collaboration, M.M.Aggarwal et.al., Phys. Lett. B403 390 (1997). 

[7] NA49 Collaboration, H. Appelshauser et al., Phys. Rev. Lett. 80 4136 (1998); C. Alt et 
al.,Phys.Rev.C68 034903 (2003). 

[8] WA98 Collaboration, M.M.Aggarwal et.al., Eur.Phys.J.C41 287 (2005). 

[9] WA98 Collaboration, M.M.Aggarwal et.al., Nucl.Phys.A 762 129 (2005) 

[10] CERES/NA45 Collaboration, J. Slivova et al., Nucl. Phys. A715 615c (2003). 

[11] STAR Collaboration, K.H. Ackermann et al., Phys. Rev. Lett 86 402 (2001); 
C. Adler et al., Phys. Rev. Lett. 87 182301 (2001). 

[12] PHENIX Collaboration, K. Adcox et al., Phys. Rev. Lett. 89 212301 (2002). S. Esumi 
et al., Nucl. Phys. A715 599c (2003). 

[13] PHOBOS Collaboration, B.B. Back et.al., Phys. Rev. Lett. 89 222301 (2002). B.Back 
et. al.,nucl-ex/0406021. 

[14] D.Molnar and S.A.Voloshin, Phys.Rev.Lett. 91 092301 (2003). 

[15] W.T. Lin et al., Nucl. Instrum. Meth. A389 415 (1997). 

[16] B. Alessandro et al., Nucl.Instrum.Meth.A493 30 (2002). 

[17] Rashmi Raniwala, Sudhir Raniwala, Yogendra P Viyogi, ALICE Internal Note; alice- 
int-2001-49. 

[18] R. Raniwala, S. Raniwala, and Y.P. Viyogi, Phys. Lett. B489 (2000) 9. 
[19] Gang Wang, Ph.D.Thesis,Kent State University, 2006. 
[20] J.Y.Ollitrault, e-print nucl-ex/9711003. 



•V 1 h 

^ 0.95 h 



0.9 - 

0.85 - 

0.8 - 
0.75 



0.75 



- £S- - -e 



-n- 



■a- 



-I. 



1 -v 

0.95 - * — 

0.9 - 

0.85 I- y hit v =0.02 
■ hit v=0.05 
0.8 h A part. v=0.02 
□ part. v=0.05 

■ ■ ■ ■ i ■ ■ ■ ■ i ■ 



-0 



-o- 



■H- 



-0 



■0- 







0.05 0.1 0.15 0.2 0.25 0.3 0.35 

X 



Figure 1: Open symbols show the ratio v c ^ jv % ^ eal for different values of x, where x = < N^us > 
I N ceil is experimentally measurable. Filled symbols show the ratio v^ ts /v l ^ eal . Squares are for v n 
= 0.05 and triangles are for v n = 0.02. The top panel is for v± and the bottom panel is for v 2 - The 
two dashed curves in each panel correspond to the different values of anisotropy v n and represent 
equation 6. The dotted curve represents equation 7. A horizontal line at the value of ratio equal 
to 1 is also drawn. 
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Figure 2: The three panels in the left column show uncorrected values of v[, RCFi and vf eom for 
charged particles and hits for an initial v 1 = 0.04 for different values of hit-occupancy x. The open 
circles are for all particles and the filled circles are for hits, as described in the text. The values are 
scaled with corresponding values for the ideal case. The dashed curve represents equation 6. The 
panels in the right column show the corresponding results for method 2, the random subdivision 
of events with the bottom column showing the value of v\ and scaled by the corresponding value 
for the ideal case. 



